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Transiting planets - lightcurve analysis for eccentric orbits 
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ABSTRACT 

Transiting planet lightcurves have historically been used predominantly for measuring 
the depth an d hence ratio of the plane t -star radii, p. Equations have previously been 
presented bv lSeager fc Mallen-Ornel as (2003) for the analysis of the total and trough 
transit lightcurve times to derive the ratio of semi-major axis to stellar radius, a/ i?*, 
in the case of circular orbits. Here, a new analytic model is proposed which operates for 
the more general case of an eccentric orbit. We aim to investigate three major effects 
our model predicts: i) the degeneracy in transit lightcurve solutions for eccentricity, 
e > ii) the asymmetry of the lightcurve and the resulting shift in the mid-transit 
time, Tmid hi) the effect of eccentricity on the ingress and egress slopes. It is also 
shown that a system with changing eccentricity and inclination may produce a long 
period transit time variation (LTTV). Furthermore, w e use our model in a reanalysis of 
HD209458b archived data bv lRichardson et a l. (2006), where we include the confirmed 
non-zero eccentricity and derive a 24/im planetary radius of Rp = 1.275Rj±0.082Rj 
(where Rj — 1 Jovian radius), which is ~ 1% larger than if we assume a circular orbit. 
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1 INTRODUCTION 

In 2000, the first transiting extr a-sol ar planet was de - 
tected bv lCharbonneau et al] I (|2000h and lHenrv et al.l (|2000l ) 
using the increasingly p rospe rous occultation method. 



ISeager fc Mallen-Ornelasl 1)2003) proposed that there are 
four principle pieces of information retrievable from any 
given lightcurve, being transit depth, AF, total transit du- 
ration, £t, transit duration between ingress and egress, 4f, 
and mid-transit time, Tmid (for which consecutive mea- 
surements may determine the orbital period). The Seager 
& Mallen-Ornelas (SMO) equations were constructed under 
the assumptions of a circular orbit and no limb darkening. 
However, as more and more eccentric transiting planets are 
discovered, the need for a more general set of equations has 
become increasingly exigent. 

From the 33 transiting planets with measured eccentric- 
ities, over 25% have eccentric orbit^J and in such cases the 
SMO equations will no longer be valid. In addition, transit 
detection probabilities favour eccentr i c syst ems by a factor 
of (1 - e 2 ) -1 , as predicted bv lBarnesI (J2007). This informa- 
tion implies that lightcurves measured for eccentric transit- 
ing planets will become common and thus a quick, efficient 
and accurate way of analyzing them is necessary. Employing 
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circularized equations to eccentric planets will clearly lead 
to unnecessary sy stematic error s of the system parameters, 
as pointed out by iBarnesI (|2007l ). 

Our intention i s to formulat e a new mo del which im- 
proves upon that by Seager & Mallen-Ornelas (2003) by al- 
lowing for orbital eccentricity. Hence, the new model must 
analytically reduce to the SMO equations in the case of 
e = 0. These equations will be able to produce sample 
lightcurves for a given set of system parameters which can 
then be compared to the real lightcurve to find a best fit. 

It will also be shown that eccentricity cannot be deter- 
mined from the primary transit lightcurve alone with cur- 
rent telescopes and any single lightcurve may be generated 
by a variety of input parameters. Another motivation for this 
work is to produce a set of analytic equations rather than 
present a computational technique for such a complex prob- 
lem. This will allow observers to produce sample lightcurves 
very quickly for any chosen system parameters and reduce 
post-observation processing time. 

We note that other authors have previo usly discussed 
the effects o f eccen tricity on transiting planets. iBurkd (|2007l ) 
and IBarnesI (J2007I ) independently focussed their discussions 
onto the probability of detection of eccentric planets, some- 
thing that we do not co nside r in this work. I n contrast, 
iTinglev fc Sackettl l|2005h and iFord et al] (|2008h presented 
models for the analysis of eccentric lightcurves. This work is 
a more elaborate development of the model developed these 
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latter authors, which we compare our model to in detail in 
§5.1. 

In §2, we briefly describe the model and the underly- 
ing assumptions (a detailed derivation may be found in the 
appendix). In §3, we discuss the implications of our model 
and degenerate solutions to a lightcurve. In §4, we show 
some typical simulations of the model for a typical sys- 
te m and in §5.1 compare it s per formance w i th th e models 
bv lTinglev fc Sackettl (|2005h and lFord et"aH (I2008D . In sec- 
tion §5.2, we discuss lightcurve asymmetry and a new effect 
where changing orbital parameters can cause a long term 
transit time variation (LTTV). In section §5.3, we employ 
our model in a rcanalysis of the HD209458b data taken by 
I Richardson et alj 1120061 ) , where we derive a new 24(im plan- 
etary radius of R P = 1.275Rj ± Q.082Rj, which is ~ 1% 
larger than assuming a circular orbit. 

In section §5.4, we discuss the effect eccentricity has on 
the ingress and egress slope shapes and conclude present- 
day telescopes cannot reliably determine eccentricity from 
lightcurve data alone. Finally, in section §6 we discuss our 
findings and conclusions. 



2 THE MODEL 

Two principle inputs for the model we present here are the 
total and trough duration times (£t and tp respectively), 
which are independent of limb darkening effects. Although, 
we stress here that the model can be coupled with the equa- 
tions and code written by Mandcl & Agol (2002) to produce 
limb darkening corrected lightcurves (see the appendix) . The 
definitions used in this paper ar e the same as that used by 
ISeager fc Mallen-Ornelaa IJ2003T ) and may be found in their 
paper, figure 1. Throughout, we make the following assump- 
tions: 

• The light comes from a single star, rather than from 
two or more blended stars. 

• The target star and planet are perfect spheres with def- 
inite edges. 

• The companion is dark. 

• The orbital period of the planet is known (from radial 
velocity surveys or consecutive Tmid times). 

• The mass of the planet (or radial velocity semi- 
amplitude) is approximately knowrQ 

Note, that we do not make the following assumptions, 
which have been adopted in previous models: 



The 



planetary 



orbit 



circular, 



ISeager fc Mallen-Orneiasl (|2003l ). 

• Mp <C M, wher e Mp is the mass of the planet, M* is 
the mass of the star. lSeager fc Mallen-Orneiasl (|2003f ). 

• The mean planet-star separation is much greater than 
both the stellar r adius, Rs, and the planetary radius, Rp, 
iFord et all (|2008l ). 

• Th e transit du r ation is much less than the orbital pe- 
riod, P. IFord et all (|2008l ). 

• The planet-star separat i on du ring the transit is near ly 
constant. Ffinglev fc Sackettl (|2005h and lFord et "ail (|2008h - 



Approximate value acceptable since M* dominates. 



A transit starts when the outer rim of the planet makes 
contact with the outer rim of the star, as seen from the 
Earth. This is defined as contact point 1 in the forthcoming 
analysis. In a similar way we define contact point 4 to be 
the exit point. We note that both points have a projected 
planet-star separation of R* + Rp (Lb), where R* & Rp 
are the stellar & planetary radii respectively. We define con- 
tact point 2 as the point when the projected planet is now 
fully inside projected star, and similar l y poin t 3 for the exit 
(figure 1 of lSeager fc Mallen-Orneiasl 1 , 200$ )). These latter 
two positions correspond to the start at end of the trough 
section of the lightcurve and have a projected planet-star 
separation of i?» — Rp (Ls). 

Our model allows for changing angular velocity as a 
function of true anomaly, as well as changing planet-star 
separation, unlike many previous analytic models. We found 
that the most accurate way to model such a system is to 
construct a full 3D geometric formulation and employ the 
conservation of angular momentum. In doing so, an elegant 
analytic expression is no longer possible and there is an argu- 
ment for a numerical approach. However, in this model, we 
do not adopt such a methodology, since analytic arguments 
offer lower computation times and the ability to reduce to 
special cases (such as e = 0, as shown in the appendix). 

In our model, we consider an elliptical orbit in a stan- 
dard frame, S, and rotate it for position of pericentre and or- 
bital inclination into a new frame, S final- This new frame 
represents what an observer from the Earth would be able 
to see, with the star centred on the origin (i.e. the projected 
orbit) . We then find the intersection points of the projected 
orbit with a circle of radius Lp and Ls respectively, which 
give us the full range of transit contact points. We are then 
able to back-transform and derive the contact point true 
anomalies, which we label as fa, fa, fa and fa. For details of 
the derivation, please see the appendix. We then convert the 
angles into times u s ing co nservation of angular momentum 
IjSertorio fc Tinettil (2002)). One significant advantage of the 
model is that we derive the projected star-planet separation, 
which c an be used as a direct input into the equations and 
code of lMandel fc Ago] (J2002J), allowing for limb darkening 
correction. Using the function D(f), which we call the dura- 
tion function, we find that the time between any two contact 
points a and b is given by: 



ta — tb 



Ha 

T 



[D(fa) ~ D(f b 



(1) 



where D(f) is defined in the appendix, J is the planet's an- 
gular momentum, fi is the reduced mass and a is the orbital 
semi-major axis 



3 IMPLICATIONS 

Any single lightcurve has three easily distinguishable char- 
acteristics; p (ratio of planet-star radii), tr and £f- Transits 
allow the measurement of orbital period, P, as well through 
the time difference between consecutive transits. For a hy- 
pothetical planet-star system, all four of these values may 
be derived using Kepler's laws and the following inputs: 



(i) Radius of the star, R, 
(ii) Radius of the planet, Rp 
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(iii) The mass of the star M* (which determines a) 
(iv) The mass of the planet, Mp (or alternatively the ra- 
dial velocity semi-amplitude, K) 
(v) The orbital inclination, i 
(vi) The eccentricity, e 
(vii) The position of pericentre, zu 

It therefore becomes clear that seven variables produce 
the resulting lightcurve shape and so there will be a degen- 
eracy in the four distinguishing features of the lightcurve, 
namely p, tr, tF and P. In other words, different combina- 
tions of system parameters can produce effectively identi- 
cal lightcurves. The implication from this is that eccentric- 
ity, e and zu cannot be determined from primary lightcurve 
data alone. Additional information is required. The model 
we present here can produce lightcurves for any given system 
parameters, but we emphasize that it is possible to produce 
ostensibly identical lightcurves from different input param- 
eters; 

ISeager fc Mallen-Ornelasl (|2003f ) appreciated that if an 
observer wishes to avoid using stellar mass-radius relation- 
ships, then one can only determine the _R* (and hence Rp) 
from a lightcurve if one has an independent and accurate 
measurement of the stellar mass. Once M* is known, _R* 
can be derived from the lightcurve data alone. We now ex- 
tend this line of thought to include e and zu. Both of these 
parameters must also be determined by an independent ob- 
servation, along with stellar mass, before the radius of the 
star can be determined. 

The ideal observation procedure would therefore be to 
use an independent measurement of M*, followed by radial 
velocity measurements to determine e, zu and Mp and then 
finally use transit observations to calculate P, R, and finally 
Rp. In the absence of an independent M, measurement, 
stellar evolution must be invoked. In section §5.3, we will 
apply a typical analysis to HD209458b, but we first compare 
our model with two other established models. 



4 VALIDATION OF THE MODEL 

For our tests of the model, we consider a Solar star, with a 
Jupiter planet at a 3-day period orbit (so a — 0.0407 AU). 
Consider that this planet to have a position of pericentre of 
20° . We now consider a range of eccentricities and produce 
the transit time durations for different inclinations. Figures 
1 & 2 show the results for multiple cases for different eccen- 
tricities, positions of pericentre and inclinations. 

Any new model proposed, must be able to be demon- 
strate that the SMO equations are a special case of the gen- 
eral form, for e = 0. For a circular orbit, we use the exact 
same assumptions as the SMO equations and therefore an 
important test of this new model is that it should reduce 
down to that form. Running a wide range of numerical cal- 
culations, we find that our model agrees precisely with that 
of the SMO equations for both tr and tF for circular or- 
bit^j. Indeed, we show in the appendix that the equations 
are analytically identical in such conditions. 
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Figure 1. Calculations for trp using our model for various e 
and i. A convergence occurs for very high e, making an accurate 
determination of i more difficult. Note how closely spaced the 
high inclinations are, making differentiation difficult here too. 
Also, we note that our model cannot generate the durations for 
orbits of e > 0.885 since the perihelion is inside the star's radius 
for such orbits. 



5 RESULTS 

5.1 Comparison & Tests of the Model 

We find our model disag rees with the model by 
iTinglev fc Sackettl (|2005h and iFord et all (|2008l ). We note 
that the models by these two authors are analytically iden- 
tical and therefore there is no need for us to compare to 
both. In figure 3, we see that there is a systematic discrep- 
ancy of around 30 seconds between Ford et al.'s model and 
that of both the model we present here and the SMO model 
for e = Qj. This systematic persists up to e « 0.5 and then 
diverges up to several minutes for very high e. We attribute 
the source of discrepancy to the different analytic expres- 
sions for the case of e = 0, where we observe that the model 
we present here reduces down to the SMO equations whereas 
that of Tingley & Sackett and Ford et al. do not, but are, 
however, simpler and quicker to utilize. 

We find that for very high eccentricities, the transit 
time duration varies very little between different inclina- 
tions, making an accurate determination of i problematic. 
In this case, we have chosen zu = 20°, but the same effect 
occurs across a range of values for w. The source of this is 
due to our test planet being in a very close orbit and so 
when we have high eccentricity, it grazes the star's surface 
and inclination becomes progressively less important. The 
effect is therefore typical but not general. 

We also observe how near zu = 270° , there is a dip in the 
transit time duration (figure 2) due to the impact parameter 
rapidly increasing as the planet enters aphelion, increasing 
the planet-star separation to a maximum; an effect which 
we expect to see. 



Note that tr and tp are independent of limb darkening 



For the case of the typical system we have described 
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Position of periastron/degrees 

e=0.2 e=0 3 e=0.4 e^0.5 — e=0.6 - 



Figure 2. Variation of tx for a typical system as a function of 
■w and e. We place a Jupiter with a 3- day orbit around a Solar 
star at i = 86° . The observed sinusoidal shape is due to the 
change in velocity at different true anomalies and the sudden dip 
at ■co ~ 270° is due to the impact parameter approaching (1 — p) 
due to increasing planet-star separation. 




Position of Periastron/degrees 

I — e=0.0 — e=0.1 e=0.2 e=0.3 — e=QA 



Figure 4. The predicted shift in Tmid for different eccentricities 
and vj, for a typical 3-day period system. The effect on Tmid 
becomes quite dramatic for very eccentric systems near 230° and 
310°. 



M a 2 [D(U) - D(h)] 
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T M i D ,true = tREF + ^-[D(f AV ) ~ D(f 1 )] (3) 
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Figure 3. Residuals in tx between the model we present here and 
Ford et al.s model, for a typical system. The difference is around 
30 seconds until we reach high e where it diverges to minutes. 



5.2 Lightcurve Asymmetry & Tmid Shifts 

Aside from solution degeneracy, the other major effect of ec- 
ce ntricity on a lig htcurve is asymmetry, as discussed recently 
bv iBarnesI (|2007l ). Since velocity is a function of phase, we 
can see that the ingress and egress slopes need not be of 
equal duration and that mid-transit time, Tmid, will be 
shifted away from the halfway point between ingress and 
egress. The direction of the shift will be a function of posi- 
tion of pericentre only, and not the sense of orbital rotation. 
The principal motivation for such a detailed investigation 
is that Tmid is predicte d to be a powerful future tool for 
exomoons dete ction (see Simon ct al. (2007)) and perturb- 
ing bodies (see lHolman fc Mu rray (2005)). We define Tmid 
shift as the difference between the apparent mid-point and 
the true mid-point. Below, we present our mathematical def- 
inition and in figure 4 we show a range of shifts for different 
system parameters in a typical system. 



ATmID = TMID,true — TMID,a 



pparent 



(4) 



where /av = 7r/2 — vj 

Transit time variation (TTV) requires Tmid to be ac- 
cura te to the second lev el for detections of perturbing bodies 
(see ISimon et al.l (J2007|)). If the target system were stable, 
then compensating for this shift would not be necessary, 
since it would only cause a constant systematic error which 
would not affect the variation of Tmid about some average. 
However, if we have a perturbing body, an eccentric tran- 
siting planet can expect to experience changing e or w over 
time. The result is that ATmid will have changed too over 
such a time-scale. 

Ergo, a perturbing body will cause two TTV sig- 
nal s, one shor t perio d signal due to the effect described 
bv lAgol et"aH (120051 ) (TTV), and one long period signal 
(LTTV) due to a combination of the expected secular 
changes and our Tmid measurement method. The difference 
between the two is subtle. The short period signal is a real 
effect, whereas the effect we describe here is a consequence 
of the measurement technique and hence a systematic error 
effect. Secular variations can cause TTV, but here the secu- 
lar variations create an additional long period signal simply 
due to the fact we typically define the mid-transit time as 
being halfway between the ingress and egress. If we always 
defined the mid-transit time as being the moment when the 
planet crosses the star's central line (Tmid, true), then this 
signal would not be present. 

In t he case of t he pr oposed, but retracted, planet 
GJ436c, iRibas et al.l IJ2008T ) suggested an important test 
would be the variation of Tmid of GJ436b, but this 
planet has a s i gnific ant eccentricity of 0.15, as measured by 
iDeming et al.l (|2007l ). For this planet, we predict ATmid is 
currently 20.95 seconds. If e, i or vo are changing over time 
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Figure 5. An example of a possible LT TV signal due to changing 
eccentricity and inclination for GJ436b. Such changes could be 
induced by a perturbing body. The effect is a result of our mea- 
surement technique plus secular variations. Over a period of a 
few years, the LTTV signal could manifest itself as transits being 
progressively late or early. 



(predicted if gravitational interactions with a third body oc- 
cur by Ribas et al.), then ATmid will also change (see figure 
5). 

In figure 4, we illustrate the shift in Tmid for the 3-day 
Jupiter test system we discussed in §4, for different e and 
w. In figure 5, we produce the long period TTV signal for 
the case of GJ436b due to the Tmid effect alone, where e is 
changing sinuoidally over a 30 year timescale with an am- 
plitude of 0.1 (values proposed by Ribas in personal com- 
munication) . We kee p position of pericentre fixed at 351° 
(jManess et al.l (|2007l )) and a l low i to change in accordance 
with that predicted bv lKozail (119621 ) to conserve angular mo- 
mentum. If zu is increasing over time, the Tmid signal would 
decrease and vice versa. Similarly, if inclination modulates 
with a larger amplitude, the Tmid signal would increase too 
and vice versa. 

As well as Tmid being shifted, the ingress and egress 
slopes can differ in duration for eccentric systems. Similar 
simulations showed that the difference between the ingress 
and egress slopes is on the order of sub-seconds to tens of 
seconds for very eccentric systems close to zu = 0° and w = 
180°. The effect is unlikely to prove useful, in general, for 
diagnosing e or zu due to typical timing errors. For example, 
in the case of the high ly eccentric system HD17156b (e = 
0.6717 as measured by iGillon et al.1 (|2007h >. the difference 
in slope times is predicted to be just 0.4 seconds. 



5.3 A Case Study for HD209458b 

Ever since its discovery, HD 209458b has been a topic 

of debate in r egard to eccentric i ty, or lack there of 

Laughlin et all i|2005l ), IWinn et al.1 (|2005l ). iKnutson et ail 



2007a|)). In real ity, no orbit can b e truly circular, and in- 
deed a study by IWinn et all (120051 ) showed that 0.0057 < 
e < .023 to within 90% confidence level. lLaughlin et al.l 
(2005) came to a similar conclusion. However, many previ- 
ous analys es of HD209458b h ave employed the SMO equa- 
tions (e.g. iTorres et al.l (J2008J)), which means the analyses 



will be somewhat flawed in assuming e = 0, when we kno w 
that this is not true, as determined by I Winn et al.l (|2005f ). 

We were interested in investigating how much difference 
including such a small eccentricity would have on the final 
result. To our knowledge, no previous authors have included 
eccentricity in the lightcurve analysis of HD209458b either, 
making such a study highly valuable. On a side note, we 
stress here a subtle point; that we do not, and indeed in 
practice cannot, use the lightcurve data alone to determine 
orbital eccentricity (see §5.4). 

For t his case analysis, we us ed the 24/j.m Spitzer re- 
sult from Richards on et al.l (| 20061 ) because the lightcurve 
is not significantly affected by limb darkening and thus 
Richardson et al. were able to derive precise timing values 
of t T = 2.979 ± 0.051 hours and t F = 2.253 ± 0.058 hours 
and a transit depth of AF = 0.01493 ± 0.00029 (see table 1 
for full list). 

Our method was to solve our equations for R* and i 
across a range of values for M*, from 0.9A/q to 1.3Mq in 
0.5A/q steps and interpolate in-between these points. The 
values and errors of the other parameters are all known and 
so in this way we may generate a lightcurve derived mass- 
radius relationship for the star (L-MRR) . We solve our equa- 
tions by employing a ne sted secant iterat ion of our equa- 
tions, using the model bv lFord et al.l ((2008) as a first guess. 
Once the optimal solution has been found, we turn our at- 
tention to calculating the error bar on the L-MRR. 

In this work, we derive our errors by employing a Monte 
Carlo bootstrap procedure with 1000 steps for selected val- 
ues of M* . This is done at each step value of A#* and then 
we once again interpolate in-between the steps to generate a 
continuous function. In each Monte Carlo run, tr, tir, p and 
e are randomly allocated according to a Gaussian probabil- 
ity distribution with the mean located at their best-known 
value and a standard deviation set to be the quoted one- 
sigma errors. The selected values and errors are shown in 
table 1. We do not consider the error in the radial velocity 
semi- amplitude, K, the orbital period, P, or the argument of 
periastron, zu, due to their very small effect on our final re- 
sult and to minimize computation time. After this stage, we 
have now produced the L-MRR and one-sigma confidence 
limits, which can be seen in figure 6. We then repeat the 
whole process but using the SMO equations (i.e. assuming 
e = 0) as a comparison. 

W e u se the eccentricit y derived by lLaughlin et al.l 
(120051 1 and lWinn et~ai] (|2005h of e = 0.014±0.005. Since the 
eccentricity is small, the effect of zu is minimal and hence 
we choose to fix it to the best known-value of zu — 83° as 

We also fix the orbital pe- 



calculated by | Winn et al 



riod, since IKnutson et al 



( 2007a]) were able to derive a very 



precise value of P = (3.52474859 ±3.8 x 10" 7 ) days. Finally, 
we fix t he radial ve locity semi-major amplitude to that de- 
rived bv lWittenmver et all (|2005l ) of K = 82.7ms" 1 , which 
is justified because both the error and the effect of the pa- 
rameter are small on our results. So the confidence limits on 
our L-MRR are derived from the errors on tr, tF, p and e. 
To find the unique solution in this case, one must em- 
ploy stellar evolution, since the errors bars on the spe ctropsi- 
cally determined mass a re quite large l|Mazeh et alj ((2000)). 
ICodv fc Sasselovl (|2002l ) provide a stellar evolution MRR by 
fitting stellar models to constant luminosity. This MRR has 
a negative gradient and so intersects our own L-MRR. We 
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Figure 6. Lightcurve derived mass-radius relation (L-MRR) for 
HD209458. The line labeled L-MRR (SMO) is the same relation 
derived assuming e = 0. The line labeled C S 02 is the stellar evo- 
lution relation derived bv \Codu & Sasselou (20031) ■ Fainter lines 
indicate one-sigma confiden ce limits. We overlay two data points, 
the lower one being that o nRichardson et al\ 1200 a) (blue) and 
the slightly higher one (red) being that of this work. 
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Figure 7. Residuals of ingress curvature of different eccentrici- 
ties compared to the e = curve. A difference of 0. 1 eccentricity 
would be typically indistinguishable with current telescopes, but 
probably not with the JWST telescope. 



find the intersection point and propagate our errors to de- 
rive our final result. The inputs and outputs of the whole 
procedure are summarized in table 1. 

We derive the 24pm planetary radius to be Rp = 
1.2757?j±0.082i?j. The result deri ved using the SMO equa - 
tions is almost identical to that of Richard son et al.l (|2006l ) 
(see table 1) and we attribute the slight differences in the 
values due to the fact we use a fixed value for K rather than 
a fixed value for M„ . 

It is found that negating eccentricity results in an sys- 
tematic underestimation of the stellar radius of HD209458 
at the 1% level, even with a very low value for e. We believe 
this outlines the importance of incorporating eccentricity 
into lightcurve analyses of other planets, even those of low 
e. 



5.4 Ingress &z Egress Slopes as a Function of 
Eccentricity 

Even in the absence of limb darkening, the ingress and egress 
slopes of the lightcurve are not straight lines but curves, 
since ^- 7^ constant during the slope part of the transit 
(where A is the occulting area). It may easily shown, that 
for any orbit, the occult ation of the star will have an S-shape 
IjMandel fc Ago! (J2002J)). We now take these equations and 
apply our new eccentric orbit model to produce the slopes. 
This will allow us to investigate the effect of eccentricity on 
such curves, which, to our knowledge, has not been investi- 
gated in detail. 

Consider the typical 3-day period system once more. If 
we set -co — 90° then the ingress and egress slope times will 
be identical and we have a symmetric lightcurve. We then 
choose i — 86° and calculate the slope shapes for different 
eccentricities. All the slopes will have the same start and 
end depth but different durations. In order to compare the 
shapes only, we scale the slope duration to go from to 1 in 
each case, in a unit we label relative time. 

Now let us consider whether it is possible to determine 
the eccentricity simply from the ingress curvature. In figure 



7, we show the residuals between the ingress curvatures for 
different eccentricities in comparison to an e = shaped 
line. We find that Ae « 0.1 gives a~5x 10 -5 change in 
normalized flux, which would be beyond the sensitivity of 
current telescopes (e.g. iBeaulieu et alj (|2008r )) but poten- 
tially feasible with the JWST missiorQ However, such an 
observation is highly unlikely to be able to better radial ve- 
locity measurements of e. 

We conclude that although it is unfeasible for current 
telescopes to distinguish between different eccentricities us- 
ing the ingress/egress slope shape or the difference between 
ingress and egress durations, future missions could do so. 
Ideally, these measurements should in the infrared, where 
limb darkening would be less significant. 



6 DISCUSSION AND CONCLUSIONS 

We have presented a new model for calculating the transit- 
ing planet lightcurves, which uses a 3D geometrical inter- 
pretation of planetary systems and may be integrated with 
the Mandel-Agol code to produce highly accurate lightcurve 
simulations. We have shown that the Seager & Mallen- 
Ornelas equations are a special case of our more general 
formulation, for when the eccentricity is zero. 

The model we present here requires two additional 
pieces of information than is needed for circular orbit 
lightcurve analysis; namely eccentricity and position of peri- 
centre. Using archived data and including the confirmed 
non-zero eccentricity of HD209458b, we derive a revised 
24pm planetary radius of 1.275±0.082i?.j , which is 1% larger 
than if we assume e = 0. We propose that the ideal method 
for determining a stellar and hence planetary radius is to 
use measurements independent of the primary transit to de- 
termine M», e, vj and Mp and then use the transits to de- 
termine P, i?* and finally Rp. In the case of no independent 
M» measurements, stellar evolution must be invoked. 

We also showed that the mid-transit time would be 
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Table 1. Summary of bootstrap inputs and results. Columns 2 & 3 are the de rived results using a bootstrap me thod and stellar evolution. 
K08 refers to the m odel presented in this pap er, SMO refers to the model by Seaaer & Malle n- Ornelak \2003i) . The final column shows 
the quoted results bv \Richardson et al\ \200a ) (R06), where we display the values as they appear in the publication. 



Parameter 



K08 



SMO 



R06 



Inputs 



P/days 


3.52474859 


3.52474859 


3.52474859 


K/ras- 1 
e 


82.7 ±1.3 
0.014 ±0.005 


82.7 ±1.3 






vj/° 


83 


- 


- 


tj- /hours 
tp /hours 


2.979 ±0.051 
2.253 ±0.058 


2.979 ±0.051 
2.253 ±0.058 


2.979 ±0.051 
2.253 ±0.058 


(Rp/R*) 


0.1222 ±0.0012 


0.1222 ±0.0012 


0.1222 ±0.0012 


Results 


M*/Mq 

R*/Rq 

P*/(kg m 3 ) 

i/° 


1 1RQ+0.096 

1 - iM -0 .079 

1 n79+ ' 055 
±.u<z_ 052 

1330 ± 230 

87.94 ±0.76 


1 174+0.096 
1 - 1 ' -0.079 

1 060+ ' 054 
l.UDU_ 051 

1390 ± 240 
88.00 ± 0.76 


1.171 
1.064 ±0.069 

88.00 ±0.85 


Mp/Mj 
Rp/Rj 


0.681 ±0.039 
1.275 ±0.082i?j 


0.690 ± 0.039 
1.261 ±0.081 


1.265 ±0.085 


p P /(kg m~ 3 ) 


407 ± 81 


426 ± 85 


- 



shifted off-centre due to lightcurve asymmetry (for eccentric 
orbits) from seconds to minutes depending upon e. Further- 
more, changing system parameters, perhaps due to a per- 
turber, could produce a long period TTV signal (LTTV) 
of the order of tens of seconds. This could manifest itself 
as progressively late or early transits over a few years time- 
scale. 

Finally, we have included an analysis on the egress and 
ingress slope curvature as a function of eccentricity and con- 
cluded that current telescopes will be unable to distinguish 
the eccentricity or position of pericentre from lightcurve data 
alone. We propose that JWST could be able to make such 
a measurement but is unlikely to match the precision pro- 
duced by radial velocity surveys. 
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APPENDIX A: DERIVATION OF THE MODEL 

Here we derive the equations underlying the model presented 
in this paper. First, consider an elliptical orbit in a frame 
S. In this frame, the Cartesian co-ordinate of the planet is 
given by: 



2 2 

x y 

h — = 1 

a 2 b 2 



(Al) 



where a = semi-major axis and b = semi-minor axis 

Now consider rotating the ellipse by an angle zu about 
the z-axis in the counter-clockwise direction. We now have 
the standard diagram for an orbit in top- view. The Si frame, 
in which the orbital inclination is 90° and the observer and 
position of periastron are defined as is historically standard. 
We also place the observer at infinite x\ position. The new 
axes are defined by: 



xi — x cos zu — y sm zu 



yi = x sin zu + y cos zu 



(A2) 
(A3) 



Now we wish to find the projected orbital view, so wc 
must change our axes a little. The standard observer would 
see the star on their LHS for zu=0, rather than the RHS as 
presently defined, so let us change the x's and y's accordingly 
and define this frame S2: 



X 2 



-Xl 



-x cos zu + y sin zu 



(A4) 



Also note: 



2/4 = J/3 + ae sin zu ■ cos i 



Z4 = Z3 + ae sm zu ■ sm 1 



(A14) 



(A15) 



By considering the points along the projected orbit 
which intersect a circle of radius L, in the £4-1/4 plane, cen- 
tred upon the stars position, we can determine the start & 
end points for the primary & secondary transits. Hence, we 
simply substitute (All) into equation (A14) and solve for 
£4: 



2,2 r 2 

x 4 + y 4 -L 



(A16) 



X4, solution = f(R,,p,i, zu,a, e) (A17) 

The resulting analytic expression for X4 is not concise 
but is possible to express with a number of substitutions, 
which we show below. Also note that L may now be replaced 
with Lb or Ls as desired. 



v = 1 — e 



(A18) 



a = 1 + - 



i(2e 2 - 4 + g) 2 + 2 cos 2 i [(2e 2 + q)(4 - q) - s\ 



(A19) 



1)2 



"J/1 



-x sm zu — y cos zu 



(A5) 



Now, let us now incline the orbit by the orbital inclina- 
tion angle i, which is done by rotating the S2 frame ellipse 
about the x-axis in a clockwise sense. This defines the S3 
coordinates as: 



X3 = X2 



J/3 = J/2 cos i + Z2 cos i — J/2 cos i 



(A6) 
(A7) 



z* 



-j/2 sini + 22 cosi 



-j/2 sini = — J/3 tani (A8) 



We substitute (A4) & (A5) into (A6) & (A7) and rear- 
range: 



-j/3 sec 1 ■ sm zu — xz cos zu 



y = X3 sm zu — j/3 sec 1 ■ cos zu 



(A9) 
(A10) 



We now substitute equations (A9) and (A10) into equa- 
tion (Al) and solve for j/3. The resulting analytic solution is 
the equation for the ellipse in transformed S3 frame. 



Q = (4? 2 )- 1 [ - aV(2 + 13e 2 ) + L 2 (7e 4 + 8v) 

+ aV cos(4i) • (8 - 6e 2 - 3q) + L 2 (e 4 - 8t>) cos(2i)l 



(A20) 



T = 



(4g 2 ) _1 [ - 8aVe 2 cos(2i)(2 + cos(2ro)) + e 2 [(-5aV 
+ 8L 2 (1 + v)) cos(2ro) + 2e 2 L 2 cos(4g7) sin 2 i\ ] 



f3 = r + g 



(A21) 



(A22) 



A = L 4 , (2avL cos if (q - 4) + 4(av cos i) 4 



-2 - 3e 2 + (2 - e 2 ) cos(2i) + 2e 2 cos(2ro) sin 2 (i) 



(A24) 



J/3a,3/3 = 

where 



i[ex 3 sin(2ro) =p \/2(l - e 2 )(qa 2 - 2x\ 



q — 2 — e + e cos(2cc) 



(AH) 



(A12) 



Now let us shift the origin to be centred about the star's 
position, the S4=Sfinal frame. 



X4 — X3 + ae cos zu 



(A13) 



e = aV[l + cos(2z)] + L 2 (g-4) 



A = — aeucos icoszu 



(A26) 



(48A) 2 • (37 2 A - 2/3-ye + 3ae 2 ) 



- 16/3(/3 2 - 9Aa) 

(A27) 



c 



4/3 2 + 12aA + 768e7g" 4 A 2 



1/3 



2-2 1/a -C 



2 2 ' s ■ (k + x/k 2 - 4C 3 ) 



(A28) 



2~ 773^1/3 



6a(« + \A 2 ~ 4C 3 )V3 

__ 64 7 2 A 2 4/3 
q 4 a 2 3a 

Finally giving: 



tin 



(A29) 
(A30) 






32 A • (32q7 3 A 2 + /37a 2 q 4 - ea 3 q 4 ) 



Ta4/2S-r±, 



a 4 9 6 VTT: 



(A31) 



The solved equation conies from solving a quartic equa- 
tion and hence we produce four solutions, which correspond 
to the start & end of the primary & secondary transits. We 
are presently only interested in the primary transit, which 
corresponds to the positive Z4 solutions. These may be dis- 
tinguished by using (A8) & (A15). Using the transform equa- 
tions, we back-transform our £4 frame solution into the S 
frame to get x(x&) and 2/(2/4). We then use simple trigonom- 
etry to acquire the true anomaly: 



iarctan ■& 
ty + arctan i? 
2-7T + arctan $ 

■where we define 



if (x — ae) > 
if (x — ae) < 
if (x — ae) > 



2/>0 
2/G5R 
2/<0 







(A32) 



(A33) 



We now take advantage of the conservation of angu- 
lar momentum and integrate (A34) between the appropriate 
limits. 



"-'t* 



(A34) 



which we integrate and then finally write the duration 
between points a and b as: 



t a — tb 






[D(fa) - D(f b )] 



(A35) 



where we define D(f) to be the duration function: 



Lightcurve analysis for eccentric orbits 9 

different times, and then input the values into the Mandel- 
Agol model, in the usual way, to produce our lightcurves. 



R* 



(A37) 



We now show that by setting e = 0, it can be easily 
shown that these equations analytically reduce down to the 
SMO equations for the timings. The substitutions reduce 
down and we find the £4 quartic solution reduces down to 
two solutions. 



XiA = — esc 1 



VL* 



a^ cos^ i 



(A38) 



i\[U- 



XiB — csc zV _L Z — a 2 cos 2 i (A39) 

Using equation (A16), we compute the 2/4 solutions too. 



2/4A 



>iy/(a 2 — L 2 ) esc 2 i 



2/4S = cos 1 



V(a 2 -^ 2 )' 



(A40) 
(A41) 



Note, that the sum of the squares of these two solutions 
is equal to L 2 , as expected. Converting these solutions back 
to the S frame we find two solutions of positive Z4: 



xa — esc i ■ \J L 2 — a? cos 2 i 


(A42) 


xb — — esc i ■ V L 2 — a 2 cos 2 i 


(A43) 


Va = Vb = -\/(a 2 - L 2 ) esc 2 i 


(A44) 



Note that once again the sum of the squares equals the 
expected result of a 2 . We therefore derive the total angular 
width of the transit to be: 



/. 



J star 



2 arcsin 



V ' L 2 — a 2 cos 2 i x 



(A45) 



Dividing this value by 27r and multiplying by the period, 
P, gives the transit durati on and we find that the result is 
the same as that derived bv lSeager fc Mallen-Ornelasl (|2003l ) 
(specifically equations (2) & (3) in the referenced paper by 
choosing L — Lb or Ls for tr and *f respectively). 

This paper has been typeset from a TfrjX/ ETfrjX file prepared 
by the author. 



D(f) = 2a/1 - e 2 arctan [ 



1 + e 



e f 

■tan 



11 



e(l — e 2 ) sin/ 
1 + e cos / 

(A36) 
We choose whether this is equivalent to Pp or tp by 
selecting L as being equal to Lb or Ls respectively, in the 
equations for X4. To employ our model with the Mandel- 
Agol limb darkening code, we simply produce the projected 
planet-star separation, z, and input it into the Mandel-Agol 
equations. In our code, we produce a table of values of z at 



